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Abstract: We reconsider in this paper the old aim of trying to understand if the observed 
realization of discrete symmetries as Parity or CP in the QCD vacuum can be satisfied from 
first principles. We show how under the appropiate assumptions implicitely done by Vafa 
and Witten in their old paper on parity realization in vector-like theories, all parity and CP 
odd operators constructed from fermion bilinears of the form ipOip should take a vanishing 
vacuum expectation value in a vector-like theory with N degenerate flavours (N > 1). 
In our analysis the Vafa- Witten theorem on the impossibility to break spontaneously the 
flavour symmetry in a vector-like theory plays a fundamental role. 
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1. Introduction 

The realization of symmetries in Nature is an essential point to understand all the phe- 
nomenology of elementary particles. Whereas space-time rotational invariance seems to be 
always realized, this is not the case for the discrete parity P, the charge conjugation, C, 
and the time reversal, T, symmetries. P and C symmetries are violated by electroweak 
processes and the only fundamental theoretical restriction comes from the CPT theorem, 
which forbids a violation of the CPT symmetry in relativistic field theory, even by the 
vacuum or ground state. 

On the other hand the large amount of experimental results on the strong interac- 
tion processes strongly suggests that P and C are preserved by this interaction, vacuum 
included. The extremely small experimental bound on the neutron dipolar electric mo- 
ment suggests that a CP violating topological term should be excluded from the QCD 
Lagrangian; and this result, combined with the solution of the U(l) problem via the axial 
anomaly has generated the well known strong CP problem [1]. Hence the understanding 
of the realization of these symmetries in QCD is one of the challenges of Theory physicists. 

One important ingredient of the fermion gauge interaction in QCD, as opposite to the 
electroweak interaction, is the fact that it is vector-like. Years ago Vafa and Witten [2] 
conjectured that parity should not be spontaneously broken in any vector-like theory as, 
for instance, QCD. 

Althought the Vafa and Witten result has been taken for a long time by the scientific 
community as a theorem, we call it conjecture for the following two reasons: 

i. V.A. and A. Galante have shown in [3] that in the Vafa-Witten work there is an 
implicit assumption: the Euclidean QCD Lagrangian in presence of an external source 
coupled to any P = — 1 local pure gluonic operator is a well defined quantum system, or at 
least a well defined statistical system, but, as a rule, this needn't be true and the system 
may not be well defined. This assumption may seem innocuous, as suggested in [4], or 
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even in [5], where the authors argued that these ill-defined theories are not physical, but 
we would like to stress that at least for a very simple statistical system, the two-dimensional 
Ising model within an imaginary external magnetic field, the theorem does not hold. Indeed 
as it was shown in [3], if parity were spontaneously broken in QCD, the theory in presence 
of an external local P = — 1 pure gluonic operator would be ill-defined. Similar concerns 
on the validity of the Vafa-Witten theorem were raised in [6], and although in [7] some 
arguments against [3] were exposed, the author claimed to have disproved the controversial 
theorem. 

%%. There are explicit counter-examples for operators built from fermion fields, as 
the parity-flavour breaking term involved in the Aoki phase [8] for Wilson fermion lattice 
QCD. Even though in the Vafa-Witten original paper it is stated in a foot note that the 
extension to fermionic operators, after the integration of fermionic degrees of freedom, is 
straightforward; this is not true: their argument can not be extended to order parameters 
constructed from fermion fields. This point, which has been stressed in [9], will also be 
discussed in this paper. 

The aim of this paper is to clarify a little bit these points and, mainly, to extend the 
Vafa-Witten conjecture to fermion bilinear order parameters, under suitable hypothesis. 

In section 2 we briefly review the Vafa-Witten argumentation, and give an alterna- 
tive way to reproduce their result, with the help of the probability distribution function 
(hereafter p.d.f.) of the order parameter. We discuss also in this section why the Vafa- 
Witten result can not be extended in a simple way to order parameters constructed from 
fermion bilinears. Since our extension of Vafa-Witten result to fermion bilinears is based 
in the p.d.f. of these operators, we give also a short review, in section 3, of the formalism 
developed in [10] to define the p.d.f. of fermionic operators. In section 4 we apply this 
formalism to the extension of the Vafa-Witten result to fermion bilinear order parameters, 
always under the assumption that a quantum theory can be consistently constructed in 
presence of these operators in the theory lagrangian. Section 5 contains our conclusions. 

2. The Vafa-Witten result: an alternative way 

The Euclidean partition function of a vector-like gauge theory can be written as the fol- 
lowing path integral over the Grassmann and gauge degrees of freedom, 



where Sym is the action for the Yang-Mills fields and the gauge-fermion coupling is a 
bilinear of the fermion fields. Ij) is the Euclidean Dirac operator and what characterizes 
vector-like theories is that the fermion determinant det^ + m), which appears in the 
effective gauge theory after integration of the fermion fields, is positive-definite for any non 
vanishing bare fermion mass. 

If we add now to the action a hermitian P = — 1 local order parameter Y(A^), which 
depends only on the gauge fields, as for instance a topological ^-vacuum term, the Euclidean 
partition function for the effective gauge theory is: 




(2.1) 
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Z = j [dA£] exp{-S YM - A J d A x Y(x) + lndet^ + m)} (2.2) 

Years ago Vafa and Witten [2] gave an argument against spontaneous P breaking in 
vector-like parity-conserving theories, as QCD. The main point in their argument was 
the observation that any arbitrary hermitian local order parameter, Y(A®), constructed 
only from Bose fields, should be proportional to an odd power of the four indices totally 
antisymmetric tensor e^ upa and therefore it would pick up a factor of i under Wick rotation. 
The local order parameter Y(A^) can therefore be written as Y(A®) = iX(A a ^) with X 
real, and the Euclidean partition function (2.2) becomes: 

Z = j [dA°] exp{-S YM ~ i\ j d 4 x X(x) + lndet^ + m)} (2.3) 

The addition of an external symmetry breaking field AY(A^) to the Lagrangian in 
Minkowski space becomes then a pure phase factor in the path integral definition of the 
partition function in Euclidean space (2.3). But a pure phase factor in the integrand of a 
partition function with positive-definite integration measure can only increase the vacuum 
energy density (free energy density). Vafa and Witten concluded that, in such a situation, 
the mean value of the order parameter should vanish in the limit of vanishing symmetry 
breaking field. 

The only objection done in [3] was that any argumentation on the vacuum energy 
density (free energy density) of the theory in presence of a parity breaking term A requires 
the previous assumption that the quantum theory is well defined in such a condition. This 
statement is less naive that what it might seem at first sight; indeed, as it was shown 
in [3], the existence of a non-parity-invariant vacuum state would imply that the theory 
in presence of a parity-breaking term is ill-defined. On the other hand, if we take this 
assumption as a requirement, the conjecture becomes a full-fledged theorem, which states 
that, if the theory with any added P = — 1 external bosonic sources is well-defined, parity 
can not be spontaneously broken. 

Let us see now an alternative way of getting the Vafa- Witten result, that makes use 
of the concept of the p.d.f. of a local operator. We do not intend to give an alternative 
proof in which we solve the issues arised in [3], but to show a different path to reach the 
same result, introducing the p.d.f. formalism. In our derivation we will assume that we 
can consistently define the quantum theory with a P-breaking local order parameter term. 

Let be Y(A^) a local operator constructed with Bose fields. The probability distribu- 
tion function of this local operator in the effective gauge theory described by the partition 
function (2.2) is 

P(c) = ^(c-i J d'x Y(x)^ (2.4) 

where V is the space-time volume and the mean values are computed over all the Yang- 
Mills configurations using the integration measure defined in (2.2). Equivalently one can 
define its Fourier transform 
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P(q) = J e tqc P(c)dc (2.5) 
which, in our case, is given by the following expression 

= ftdAft eMSYM + (g ~ A) J d 4 x r (x) + In detff + m)} 

W /[d^]exp{-5y M - A/d 4 x F(x) + lndet$ + m)} l ' j 

or, in a short notation: 

P(q) = (eMy J d 4 x Y(x))^ (2.7) 

If we calculate P(c) for a local P-breaking operator in absence of a P-breaking term 
in the action of (2.2) (A = 0) we expect a Dirac S distribution centered at the origin if the 
vacuum state is non degenerate. On the contrary, if the vacuum state is degenerate, the 
expected form for P(c) is: 

P(c)=^2w a 5(c-c a ) (2.8) 

a 

where c a is the mean value of the local order parameter in the vacuum state a and w a are 
positive real numbers which give us the probability of each vacuum state (J2 a w a = !)• 

If the only reason to have a degenerate vacuum is the spontaneous breaking of a discrete 
Z2 symmetry, we expect in the more standard case two symmetric vacuum states (±c a ) 
with the same weights w a ( in the most general case an even number of vacuum states with 
opposite values of c a ). The probability distribution function P(c) will be then the sum of 
two symmetric Dirac 5 with equal weights: 



P(c) = ^5(c-c a ) + ^5(c + c a ) (2.9) 



and its Fourier transform 



P{q) = cos(qc a ) (2.10) 

Now we are ready to apply this formalism to the case we are interested in. The relevant 
fact now, as in the Vafa-Witten paper, is the fact that the local P-breaking order parameter 
is a pure imaginary number ^(^4^) = iX(A^)). In such a case we get for P(q): 

P(q) = (eM~ I d A x X{x))} (2.11) 

but the mean value of a real and positive quantity computed with a well defined proba- 
bility distribution function (positive integration measure) can not be negative. Were the 
symmetry spontaneously broken we should get for P(q) a cosine function (2.10) (a sum 
of cosines in the most general case) which takes positive and negative values. Since this 
case is excluded, P(q) should be a constant function equal to 1, representing a symmetric 
vacuum state. 
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The generalization of the Vafa-Witten argumentation to fermionic bilinears upon in- 
tegrating out the fermion fields, as suggested in [2], can not be carried out, as discussed in 
[9]. Indeed if we add a fermionic bilinear P = — 1 term, Xitp^tp, to the fermionic action, 
we get in the effective Yang-Mills theory an extra-term to the pure gauge action with a 
non-vanishing A-dependent real part. Thus the contribution to the effective gauge theory 
of the symmetry breaking term in the integrand of the partition function is no longer a 
pure phase factor and therefore how this term modifies the vacuum energy-density cannot 
be stated a priori. In fact there is a well known case, the Aoki phase [8] in lattice QCD 
with two degenerate flavours of Wilson fermions, in which notwithstanding the integration 
measure is non negative, Parity is spontaneously broken. 

3. The p.d.f. for fermionic local operators 

Our aim is to extend, as much as possible, the Vafa-Witten result for pure gluonic operators 
to fermion bilinear local operators. To this end we will make use, as in the previous section, 
of the p.d.f. of local operators constructed with Grassmann fields. Since this formalism is 
not, in our opinion, very standard, we want to devote this section to summarize the main 
results reported in [10] on this subject. 

The motivation to develop this formalism was precisely the study of the vacuum in- 
variance (non invariance), in quantum theories regularized on a space time lattice, under 
symmetry transformations which, as chiral, flavour or baryon symmetries, involve fermion 
fields. The numerical determination of the p.d.f. of the order parameter is a standard 
procedure when analyzing spontaneous symmetry breaking in spin systems or in quantum 
field theories with Bose fields. Indeed these degrees of freedom can be simulated directly 
on a computer. However in the numerical simulation of a Lattice Gauge Theory (LGT) 
with dynamical fermions the Grassmann fields, which at present can not be simulated in a 
computer, must be integrated analytically. Then even if the ground state of the theory is 
not invariant under the chiral, flavour or baryon symmetries, in the analytical procedure 
we integrate over all possible vacuum states, the order parameter vanishing always inde- 
pendently of the symmetry realization. This is the reason why the standard procedure to 
study, for instance, the chiral properties of the QCD vacuum in numerical simulations on 
a lattice is to add a mass term to the Euclidean action, which breaks chiral symmetry and 
to compute the chiral condensate {ipip) as a function of the mass m performing, at the end, 
a numerical extrapolation to the chiral limit in order to see if takes a non-vanishing 
value in this limit. 

Notwithstanding that Grassmann variables cannot be simulated in a computer, it was 
shown in [10] that an analysis of spontaneous symmetry breaking based on the use of the 
p.d.f. of fermion local operators, and therefore free from extrapolations in the symmetry 
breaking field, can also be done in QFT with fermion degrees of freedom. 

The procedure is similar to the one used in the previous section for bosonic degrees 
of freedom. The starting point is to choose an order parameter for the desired symme- 
try 0(ifj,ip) (typically a fermion bilinear) and characterize each vacuum state a by the 
expectation value c a of the order parameter in the a state 
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c a = i J (0(x)) a d 4 x (3.1) 
The p.d.f. P(c) of the order parameter will be given by 

P(c) =Y^w ol 5{c-c a ) (3.2) 

a 

which can also be written as [10] 

P(c) = (s(y J 0{x)d\-c)^ (3.3) 

the mean value computed with the integration measure of the path integral formulation of 
the Quantum Theory. 

The Fourier transform P(q) = J e iqc P(c)dc can be written, for the theory described 
by partition function (2.1), as 

p, JidA^d^expi-SYM + ^ + m^ + ^Jd^x O(x)} 

In the particular case in which O is a fermion bilinear of ip and ■0 

0(x) = ${x)d$(x) (3.5) 
with O any matrix with Dirac, color and flavour indices, equation (3.4) becomes 

= f[dA^][dHiP] exp{-S YM + $<p + m + f6)^} 
[<1) J[dA*][d$d4]exp{-S YM + $P + mW} [ '> 
Integrating out the fermion fields in (3.6) one gets 

= /[^le-^det^ + m + fO) 
[q) f[dAfa- s YM det$ + m) { ' 

which can expressed also as the following mean value 

computed in the effective gauge theory with the integration measure 



[cM^e -5 ™ det$ + 



m) 



The particular form expected for the p.d.f. P(c), P(q), depends on the realization of 
the corresponding symmetry in the vacuum. A symmetric vacuum will give 

P(c) = 5(c) (3.9) 
P(q) = 1 
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whereas, if we have for instance a £7(1) symmetry which is spontaneously broken, the 
expected values for P(c) and P(q) are [10] 




(3.10) 



the last being the well known zeroth order Bessel function of the first kind, Jo(qco). 

In the simpler case in which a discrete Z{2) symmetry, as Parity, is spontaneously 
broken, the expected form is 



or a sum of symmetric delta functions (-P(c)) and a sum of cosines (P(q)) if there is an 
extra vacuum degenerancy. 

This formalism has been successfully applied to the analysis of spontaneous chiral 
symmetry breaking in LGT with Kogut-Susskind fermions [10], as well as to find a new 
phase in finite baryon density two-colors QCD with diquark condensation and spontaneous 
breaking of baryon symmetry [11]. 

The analysis of the diquark condensation phase in finite density two-color QCD with 
four Kogut-Susskind flavours done in [11] showed also one of the potentialities of this 
approach. Indeed the order parameter used in the search for diquark condensation phase 



when added to the SU(2) Lagrangian as an external source that breaks baryon symmetry, 
place us against a tedious problem, the well known sign problem. After integrating out 
the fermion fields we get the Pfaffian of a matrix and this Pfaffian is not positive definite. 
The standard approach to perform numerical analysis of spontaneous symmetry breaking, 
based on the computation of the condensate at non- vanishing external source, does not work 
because standard Monte Carlo techniques apply to systems with a well defined (positive 
definite) Boltzmann factor. The p.d.f. of the diquark condensate is computed, on the 
contrary, at vanishing external source, and the sign problem is absent in this case. 

We are ready now to apply this approach in the next section to the analysis of spon- 
taneous P breaking in vector-like theories, by investigating order parameters constructed 
from bilinears of the fermion field. 

4. Vacuum expectation values of P — — 1 fermion bilinears 

In this section we will show that all local bilinear P = — 1 gauge invariant operators of the 
form t/jOip, with O a constant matrix with Dirac, color and flavour indices, should take a 
vanishing vacuum expectation value in any vector-like theory with N degenerate flavours. 



P( c ) = ^5{c-c ) + ^5{c + c ) 



P(q) = cos(gco) 



(3.11) 
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Let us start with the one-flavour case since, as it will be shown, it is a special case. The 
standard hermitian local and gauge invariant P order parameter bilinear in the fermion 
fields is: 



tyOi\) = iip'js'ip 

Now we can apply the result of the previous section, and in particular equation (3.8), 
which give us the p.d.f of O in momentum space 

_ /detP + m+faU 
W \ det^ + m) / V 7 

The determinant of the Dirac operator in the denominator of (4.1) is positive definite, 
but the numerator of this expression, even if real, has not well defined sign. The final 
form for P(q) will depend crucially on the distribution of the real eigenvalues of 75^ + m). 
Therefore we cannot say a priori whether P(q) will be the constant function P(q) = 1 
(symmetric vacuum) or a cosine function (spontaneously broken P). 

Let us go now to the N flavour case {N > 1). The most general P = — 1 hermitian 
and gauge invariant local order parameters ipOi/) which can be constructed are: 



ivp'jsTtp (4.2) 

with r any of the hermitian generators of the SU(N) flavour group. However, since flavour 
symmetry cannot be spontaneously broken in a vector-like theory [12], we will restrict our 
analysis to the flavour singlet case. 



#75^ = i^ulbTpu + i^dlb^d + i^sl5^s + ••• (4.3) 

Let us assume that {i'tpulhi'u) = ±co / 0. Since flavour symmetry is not spontaneously 
broken we should have: 



{i^ul^u) = {i^dlb^d) = (#s75^s) = ••• (4.4) 

Thus the system will show two degenerate vacua with all the condensates oriented in 
the up (down) direction. 

We want to remark that the "ferromagnetic" nature of these two vacua is imposed 
by the realization of flavour symmetry in the vacuum. Otherwise one could imagine also 
"antiferromagnetic" vacua with antiparallel condensates, or even more complex structures. 

In vector-like gauge theories the interaction between different flavour is mediated by 
gluons. Then, assuming {i^ulb^u) 7^ the impossibility to break flavour symmetry [12] 
could suggest that the gauge interaction would favor "ferromagnetic" vacua with parallel 
oriented condensates. However the actual dynamics can also be more complex. Indeed 
we will argue now that a non vanishing condensate (iipujsipu) 7^ 0, which would imply 
spontaneous breaking of P, CP, T and CT, can be excluded. 
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Let us apply equation (3.8) to the computation of the p.d.f. in momentum space P u d(q) 
of i^ulb^u + iipdlbipd- 



Pud{q) = 




det^ + m-f 75 ) \ 
det$ + m) ) 



(4.5) 



where-0 + m in (4.5) is the one flavour Dirac operator and the mean value is computed in 
the theory with N-degenerate flavours. 

Equation (4.5) give us the mean value, computed with a positive definite integration 
measure, of a real non- negative quantity. Thus P u d(q) will be a positive definite, or at 
least a non negative definite, quantity. But were {i^ulb^u) / Owe should expect a cosine 
function for P u d(q) since i^ulb^u and i^dlb^d are enforced to take the same v.e.v. for 
flavour symmetry. Since the positivity of P u d(q) excludes such a possibility, we can conclude 
that all the pseudo-scalar condensates iV'/TsV'/ take vanishing expectation values. 

In order to give a simple illustration of some of this ideas, one can think of a physical 
system composed by two replicas of the two-dimensional Ising model. The spins s a of 
replica a and s b of replica b live on the sites of a two-dimensional lattice. They are coupled 
inside each replica with the standard nearest-neighbor ferromagnetic coupling, and we add 
an ultra-local interaction between s a and s b at the same site. 

The Hamiltonian of this system can be written in the following way: 



This system has the standard Z{2) P symmetry ( H is invariant under the change of 
sign of all spins) plus an extra Z{2) symmetry, the replica symmetry, as a consequence of 
the invariance of TL under the permutations of the a and b indices. 

Let us neglect for a while the replica-replica interaction (k = in (4.6)) and choose 
F = jl?p > 0.44 (ferromagnetic phase). If we denote the magnetization per site of replica 
a and b as m a and m&, we can construct two order parameters: (m a — nib) and (m a + nib). 
Both are P = — 1 order parameters and the first one is also order parameter for the Z(2) 
replica symmetry. 

Under the previous assumed conditions {F > 0.44, k = 0) we have four degenerate 
vacua (++, , H — , — h) characterized by the following values of the order parameters: 



If we switch on now the replica-replica interaction (A; ^ 0) we expect that the vacuum 
degeneracy will be reduced. Indeed k > values will enforce a parallel orientation of spins 
of different replicas, whereas k < will favor an anti-parallel orientation. In the first case 
we will get two "ferromagnetic" degenerate vacua preserving the replica symmetry. In the 




(4.6) 



<i,j> 



(m a + m b ) ++ = 2m (m a - m b } ++ = 
{m a + m b ) — = -2m (m a - m b ) — =0 

(m a + m b )+- = (m a - m b ) + - = 2m 
(m a + m 5 }_ + = (m a - = -2m 
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second case we get two "antiferromagnetic" degenerate vacua with spontaneous parity and 
replica breaking. 

What we can learn from this simple model is that a small replica-replica interaction is 
enough to break some of the vacuum degeneracy. Then "ferromagnetic" or "antiferromag- 
netic" degenerate vacua are selected, depending on the ferro- antiferromagnetic character 
of the replica-replica interaction. 

5. Conclusions 

We have reconsidered in this paper the old aim of trying to understand if the observed 
realization of discrete symmetries as Parity or CP in the QCD vacuum can be satisfied 
from first principles. Although this aim is at the moment too ambitious, we have shown 
how under the appropiate plausible assumptions implicitely done in [2], all parity and CP 
odd operators constructed from fermion bilinears of the form -ipOijj should take a vanishing 
vacuum expectation value in a vector-like theory with N degenerate flavours (N > 1). 
In our analysis the Vafa-Witten theorem on the impossibility to break spontaneously the 
flavour symmetry in a vector-like theory [12] plays a fundamental role. Indeed our re- 
sult does not apply to lattice QCD with two degenerate Wilson fermions, where a phase 
with spontaneous Parity breaking has been found [8]. It follows because, notwithstanding 
the integration measure is non negative, the Vafa-Witten theorem on the impossibility of 
breaking spontaneously a vector symmetry in a vector-like theory is not applicable to this 
case, for the Wilson quark propagator is not bounded in an arbitrary background gauge 
field. Indeed flavour symmetry is spontaneously broken in the Aoki phase. 
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